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ABSTRACT 

We report on the results of a 50 ks BeppoSAX observation of 1RXS J 170849-4009 10, one of the five (plus a 
candidate) known anomalous X-ray pulsars. The BeppoSAX data are consistent with a power-law plus black- 
body spectral decomposition, making 1RXS J 170849-4009 10 the fourth source of this class for which such a 
spectral decomposition was found. The inferred power-law slope and blackbody temperature are V ^2.6 and 
kTfiB ~0.46 keV, respectively. We found significant energy-dependence of the pulse profile, a remarkable feature 
for an AXP. By using the power-law plus blackbody decomposition we detected a significant variation in at least 
one spectral parameter, the power-law photon index, as a function of the pulse phase. This is the first significant 
detection of spectral parameter variation in an AXP. The implications of these results are briefly discussed. 

Subject headings: stars: neutron — pulsar: individual (1RXS J 170849-4009 10) — pulsars: general — X-ray: 
stars 



1. INTRODUCTION 

After more than 20 years since the discovery of pulsations 
from IE 2259+586 the nature of the Anomalous X-ray Pulsars 
(AXPs) is still an open issue. Although we can be reasonably 
confident that AXPs are magnetic rotating neutron stars (NSs), 
their energy production mechanism is still uncertain. It is also 
unclear whether they are solitary objects or are in binary sys- 
tems with very low mass companions (for a review see Israel 
et al. 2001 and references therein). Different production mech- 
anisms for the observed X-ray emission have been proposed, 
involving either accretion or the dissipation of magnetic energy. 
Among the properties that distinguish AXPs from known mag- 
netic (> 10 12 G) accreting X-ray pulsars found in High and 
Low Mass X-Ray Binaries (HMXBs and LMXBs) are: (i) spin 
periods in a narrow range (^6-12 s), (ii) very soft and absorbed 
X-ray spectra, (iii) relatively stable spin period evolution, with 
long term spin-down trend, (iv) flat distribution in the Galac- 
tic plane and three clear associations with supernova remnants. 
There are currently five ascertained members of the AXP class 
plus one likely candidate. 

The relatively bright source 1 RXS J 1 70849-4009 1 was dis- 
covered early in the ROSAT mission (Voges et al. 1996); how- 
ever only in 1997 this source attracted much attention because 
of the ASCA discovery of ~ 1 1 s pulsations (Sugizaki et al. 
1997). Based on the pulse period and unusually soft X-ray 
spectrum the source was tentatively classified as a candidate 
AXP. This interpretation was confirmed through ROSAT High 
Resolution Imager (HRI) observations which provided the first 
measurement of the period derivative P^2xlO -11 ss" 1 (Israel 
et al. 1999a). 1RXS J 170849-4009 10 is one of the two AXPs 
for which a phase-coherent timing solution was obtained by a 
systematic monitoring program with the RossiXTE PC A (Kaspi 
et al. 1999). The source was found to be a quite stable rotator 
with phase residuals of only ~1%, i.e. comparable to or smaller 
than those measured for most radio pulsars. However, on 1999 
September the RossiXTE satellite detected a sudden spin-up 



event from 1RXS J 170849-4009 10 which was interpreted as a 
"glitch" similar to those observed in the Vela and other young 
radio pulsars (Kaspi et al. 2000). 

A search for optical counterparts in the field of 
1RXS J 170849-4009 10 was carried out by Israel et al. (1999a). 
These authors, based on two refined ROSAT HRI positions, 
were able to rule out a massive early type companion (a distant 
and/or absorbed OB star would appear more reddened). How- 
ever, the images were taken from a 1.5 m telescope and were 
not deep enough to constrain any other proposed theoretical 
scenario such as a low-mass companion, a residual disk or a 
magnetar. 

An association between 1 RXS J 170849-4009 10 and the su- 
pernova remnant (SNR) G346.6-0.2 located ^12' away was 
proposed by Marsden et al. (2001). However, as discussed 
by Gaensler et al. (2001) such an association appears to be 
unlikely. An image at 1.4 GHz showed the presence of an 
arc of diffuse emission -8' away from 1 RXS J 170849-4009 10 
(Gaensler et al. 2001), which was interpreted as a previously 
unknown supernova remnant (G346.5-0.1). Also in this case 
there are no convincing arguments for a physical association 
between G346.5-0.1 and 1 RXS J 170849-4009 10. No radio 
emission was detected from 1 RXS J 170849-4009 10, with up- 
per limits of 3 mJy on the continuum (5er at 1.4 GHz; Gaensler 
et al. 2001). In this paper we report the results obtained from 
a BeppoSAX observation of 1RXS J 170849-4009 10 that took 
place before the "glitch"-like event detected by RossiXTE. A 
two-component spectrum, i.e. a power-law plus a black body, 
was found. Moreover the BeppoSAX observation revealed 
significant energy-dependent pulse profile. Pulse phase spec- 
troscopy shows that, in the context of the blackbody plus power 
law spectral decomposition, this effect is likely connected to 
a photon index variation. This observation therefore provided 
the first evidence for pulse phase-related spectral variations in 
an AXP. Implications of these results are also briefly discussed. 
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FIG. 1.— 1RXS J 170849-4009 10 MECS and LECS light curves (left pan- 
els labelled from (a) to (d)) folded to the best period (P=10.99915 s) for four 
different energy intervals. For clarity two pulse cycles are shown. Zero phase 
was (arbitrarily) chosen to correspond to the minimum in the 0.3-2 keV folded 
light curve. The results of the pulse phase spectroscopy are shown for selected 
free parameters (right panels; the absorption column in units of 10 22 ). Phase 
intervals refer to those shown in the left panel. 

2. OBSERVATIONS AND RESULTS 

BeppoSAX observed 1RXS J 170849-4009 10 between 1999 
March 3 1 14:27 UT and April 1 22:50 UT with imaging Narrow 
field Instruments: the Low-Energy Concentrator Spectrometer 
(LECS; 0.1-10 keV; Parmar et al. 1997; 26 ks effective expo- 
sure time) and the Medium-Energy Concentrator Spectrometer 
(MECS; 1.3-10 keV; Boella et al. 1997; 52 ks effective expo- 
sure time). 

2. 1 . Timing analysis 

The arrival times of the 0.1-10 keV photons from 
1RXS J 170849-4009 10 were corrected to the barycenter of the 
solar system and a 1 s binned light curve was accumulated. The 
MECS counts were used to accurately determine the pulse pe- 
riod. The data were divided into 12 time intervals, and for each 
interval the relative phase of the pulsations was determined. 
These phases were then fit with a linear function giving a best- 
fit period of 1 1.999 15 ±0.00002 s. The background subtracted 
light curves, folded at the best period in different energy ranges 
(Figure 1 ; left panels) show an energy-dependent profile. In 
particular the phase interval of the minimum in the energy in- 
terval 0.3-2keV folded light curve corresponds to a maximum 
in the 6-10keV folded light curve. Moreover the pulsed frac- 
tion (semi-amplitude of modulation divided by the mean source 
count rate) decreases from ^40% to ^30% from the lowest to 
the highest energy band. 

2.2. Spectral analysis 

PHA spectra were obtained from the BeppoSAX position 
(RA. = 17 h 08 m 48 s , Dec. = -40°08'58"; equinox 2000; 90% 
confidence level radius of 22") of 1RXS J 170849-4009 10 us- 
ing an extraction radius of 4' and 8' for the MECS and LECS, 
respectively. Background subtraction was performed using 
both standard blank field exposures and background regions 
taken from the observation of 1RXS J 170849-4009 10 far from 
the position of the AXP ending up with similar results. The 
PHA spectra were rebinned in order to have more than 40 



counts in each bin such that minimum x 2 fitting techniques 
could be reliably used. All those bins which were consistent 
with containing zero counts after background subtraction were 

Table 1 

BeppoSAX PHASE AVERAGED FIT OF 1RXS J170849-400910. 



Spectral Parameter 


PL 


PL + BB 


BB + BB 


N H (10 22 atom cm" 2 ) 


1.88±0.08 


1.42±0.15 


0.9±0.1 


r 


3.28±0.05 


2.62±0.17 




PL flux 


4.26 


3.30 




kT m (keV) 




0.46±0.03 


0.50±0.02 








1.54±0.06 


BB radius (km) 




4.0±0.4 


4.4±0.2 








0.30±0.02 


BB flux 




1.2 


2.8+1.3 








1.7 


X 2 /dof 


1.24 


1.00 


1.10 


Lx (10 35 erg s~') 


8.5 


3.6 


2.0 



NOTES — Fluxes and luminosities refer to the 0.5-10 keV band (x 
10~" erg s~'cnT 2 . Fluxes are not corrected for the interstellar ab- 
sorption. Flux uncertainties are about 10%. The source luminosities 
were derived by setting N H = and assuming a distance of 5 kpc. 



rejected. Moreover the analysis of the MECS and LECS spectra 
was restricted to the 1.8-10keV and 0.1-5keV ranges, respec- 
tively, where the calibration of response files is more accurate. 
A constant factor free to vary within a predetermined range was 
applied in the fitting to account for known normalization differ- 
ences between LECS and MECS. 

A simple power-law model did not fit the data well (reduced 
X 2 of 1.24 for 204 degrees of freedom, hereafter dof). All other 
single component models that we tried produced even worse 
results. A much better fit (see Figure 2) was obtained by in- 
cluding a soft thermal component, a blackbody, in addition to 
the power-law (reduced x 2 of 1.00 for 202 dof). This model 
was successfully fitted to the spectra of three other AXPs. An 
F-test shows that the inclusion of the blackbody component is 
highly significant (probability of ^6cr). The best fit was ob- 
tained for an absorbed (N H = (1.42±0.15)x 10 22 atom cm" 2 ) 
power-law with photon index T = 2.6±0.2 and a blackbody 
component with temperature of kT BB = 0.46±0.03keV (90% 
c.l. reported; see Table 1). The unabsorbed 0.5-10keV flux 
was 1.3 x 10" 10 erg s" 1 cm" 2 . The blackbody component ac- 
counts for about 36% of the total absorbed flux in the same 
band. Figure 2 shows the spectral shape and components of 
1RXS J 170849-4009 10 as determined by BeppoSAX. We also 
tried to fit the spectra using different spectral decompositions. 
Among these we find a relatively good fit (reduced x 2 of 1.10 
for 202 dof) with a two blackbody model (see Table 1). It is 
worth mentioning that also in this case we obtained a soft com- 
ponent (i.e. a black body) with a characteristic temperature sim- 
ilar to that inferred with the power-law plus black body model. 

The data from the High Pressure Gas Scintillation 
Proportional Counter (HPGSPC) and the Phoswich De- 
tector System (PDS) did not provide any useful in- 
formation on 1RXS J 170849^009 10. In fact, due 
to the large FOVs of these instruments, the relatively 
short exposure time (~25ks), and steep spectrum of 
1RXS J 170849-4009 10 the source was not detected. 
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FIG. 2.— LECS and MECS energy spectra of 1RXS J 170849-4009 10. The 
residuals (in units of \ 2 ) of the best fit are also shown (see the text for details). 
The power-law and blackbody components are shown with dotted-dashed and 
dashed lines, respectively. 

2.3. Phase resolved spectroscopy and pulse profiles 

Pulse phase spectroscopy (PPS) was carried out with the 
MECS and LECS data. A set of five phase-resolved spec- 
tra (phase boundaries 0.0, 0.2, 0.4, 0.6, 0.8) were accumu- 
lated. After rebinning and background subtraction these were 
then fit with the power-law plus blackbody model described 
in Sect. 3.1. Due to the small number of photons in the 
LECS spectra, we removed all the PHA channels below 1 keV. 
Initially, the blackbody temperature was fixed at its phase- 
averaged best-fit value and only the power-law parameters and 
blackbody normalisation were allowed to vary giving a cumu- 
lative (over the whole set of spectra) reduced y} of 1.10 for 
466 dof. The fits were then repeated with the power-law com- 
ponent fixed and blackbody parameters free to vary resulting 
in a reduced x 2 of 1.32 for 466 dof for the best fit. Then all 
the parameters were varied and fitted together. In the latter 
case the best fit gave a reduced \ 2 of 1.04 for 456 dof. An 
F-test shows that the freeing of the blackbody and power-law 
phase-averaged parameters is highly significant (the reduction 
in x 2 has a formal probability of 9x 10" 4 and 1 x 10" 19 for the 
blackbody and power-law parameters, respectively). In Fig- 
ure 1 (right panels) and 3 the results of the PPS are shown for 
the most interesting parameters and phase intervals. 

A check was performed by fitting the power-law photon in- 
dex values obtained in each phase intervals with a constant (see 
right panels of Figure 1). This set a ^3c significant variation in 
the power-law photon index parameter. Statistical uncertainties 
prevent a firm detection of variations in the other parameters 
(see right panels of Figure 1). 

The upper two panels of Fig. 4 show the pulsed fraction 
versus energy during BeppoSAX observation for the first two 
harmonics. These values were obtained by fitting the corre- 
sponding light curves with two sinusoidal functions. Note that 
definitions of the pulsed fraction which involve the maximum 
and the minimum of the folded light curve should be used with 
care as they are dependent on the binning time and, therefore, 
the presence of unusually low or high data points. The sinu- 
soidal fit, when feasible, is less sensitive to these data points. 
The first harmonic decreases from ^36% to ^26% as the en- 



ergy increases from 0.5-2 keV to 6-10keV. A constant value of 
~10% is inferred for the second harmonic. The last two panels 
refer to the ratio between the power-law and the total absorbed 
fluxes (i.e. the blackbody plus power-law spectral model; third 
panel), and the ratio between the blackbody and the power- 
law unabsorbed fluxes (lowest panel). From the comparison of 
these quantities we can infer that: (i) there is evidence for a 
decrease of the 1" harmonic pulsed fraction at energies above 
5keV, although the statistics are poor; (ii) there is evidence 
for an anti-correlation between the power-law component and 
the 1" harmonic pulsed fraction (see first and third panels), 
and (iii) evidence for a direct correlation between the latter 
and the blackbody component (see first and fourth panels). 
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FIG. 3.— LECS and MECS spectra of 1RXS J 170849-4009 10 for two se- 
lected phase intervals (0.2-0.4) and (0.8-1.0) as reported in Figure 2. The 
power-law photon index variation is clearly evident. 

3. DISCUSSION 

The BeppoSAX spectrum of 1RXS J 170849-4009 10 is 
well modeled by the sum of a (absorbed) relatively 
steep power-law and a low-energy blackbody. Therefore 
1RXS J 170849-4009 10 is the fourth AXP, after 4U0142+614 
(White et al. 1996; Israel et al. 1999b), IE 2259+586 (Cor- 
bet et al. 1995; Parmar et al. 1998) and IE 1048.1-5937 
(Oosterbroek et al. 1998), for which such a spectral decom- 
position has been detected. We note that these sources are 
also the ones for which good spectral data are available. This 
possible decomposition for 1RXS J 170849-4009 10 was first 
suggested by Sugizaki (1997) although not statistically sig- 
nificant using the ASCA data. The comparison between the 
BeppoSAX and ASCA observations reveal that the source has 
remained nearly at the same (absorbed) flux level (~4.4 x 10~ n 
erg cm" 2 s" 1 and ~4.3xl0" n erg cm" 2 s" 1 in the 0.8-10keV 
band for BeppoSAX and ASCA, respectively). Also the spec- 
tral parameters (in the blackbody plus power-law model) are 
similar to those found by ASCA with the exception of the tem- 
perature of the blackbody which is higher in BeppoSAX: this is 
not unusual as BeppoSAX instruments have a higher sensitivity 
below 1 keV allowing a better evaluation of the absorption col- 
umn and, therefore, of the blackbody component parameters. 

We also tried to fit the spectrum of 1RXS J 170849-4009 10 
with other models. Among these we used two blackbodies. 
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FIG. 4. — Pulsed fraction of the first two harmonics of 
1RXS J 170849-4009 10 as a function of the energy (upper two panels), 
together with the power-law to total absorbed flux ratio and the blackbody to 
power-low unabsorbed flux ratio (lower two panels). See text for details. 

Such a decomposition is in agreement with accreting, magnetic 
field decay and cooling models which have been invoked to ac- 
count for the main physical mechanism(s) responsible for the 
X-ray emission of AXPs (see also Thompson & Duncan 1996 
and Heyl & Hernquist 1997). Two (or more) concentric regions 
with different temperatures are in fact expected: the innermost 
corresponding to the polar caps of the neutron star (which is 
also the hottest) and a larger area around the magnetic caps 
characterised by a lower temperature (see also DeDeo et al. 
2000; Ozel et al. 2001; Perna et al. 2001). Regardless the 
origin of the emission from these regions, we note that the size 
of the blackbodies is always smaller than the neutron star sur- 
face, even for an unrealistic distance to the AXP of ~45kpc, 
with the smallest region (R B b ~0.3 km) being also the hottest 
(kT ~1.5 keV). The results of a more detailed and systematic 
spectral and timing study of the AXPs observed by BeppoSAX 
will be reported elsewhere. 

An analysis of the pulse shape of 1RXS J 170849-4009 10 
as a function of energy (in the 0.3-10keV) reveals a promi- 
nent variation with the minimum at low energies correspond- 
ing to the maximum at high energies (see also Figure 1). Cor- 
respondingly, pulse phase spectroscopy detected a significant 
variation for at least one spectral component, the power-law 
photon index, as a function of phase. The peak in the pulse 
profile at highest energy also corresponds to the lowest value 
of the photon index, similar to what observed in accreting X- 



ray pulsars (Makishima et al. 1999 and references therein). 
Although energy-dependent changes in the pulse shape were 
already observed in the past for 4U0142+614 (White et al. 
1996; Israel et al. 1999b; Paul et al. 2000), the variations of 
1RXS J 170849-4009 10 are also accompanied by a nearly total 
phase reversal between the low and high energy band, and are 
likely due by a changing power-law slope (as suggested by Su- 
gizaki et al. 1997). We note that, in the past, the lack of any 
conspicuous change in the pulse profiles as a function of en- 
ergy and/or pulse phase-resolved spectra of AXPs was used to 
argue against the possibility that these sources are accreting X- 
ray objects. We finally note that the pulse fraction as a function 
of energy shown in Figure 4 is not in disagreement with that 
recently reported by Ozel et al. (2001) which used a different 
definition of the fractional contribution (and assuming only one 
harmonic) to the pulsations. 

1RXS J 170849-4009 10 is so far the only AXP for which a 
sudden spin-up was observed (Kaspi et al. 2000). This was 
interpreted as a glitch similar to that observed in the Vela and 
other young radio pulsars. However, we note that glitches could 
in principle be detected also in accreting (spinning-down) X- 
ray sources with a sufficiently high magnetic field strength 
(in analogy with the known characteristics of radio pulsars) 
if they are in a low noise level phase, as indeed AXPs are 
known to be. A way to distinguish, in the near future, whether 
1RXS J 170849-4009 10 experienced a radio pulsar-like glitch 
or, perhaps, an accreting X-ray pulsar-like spin-up behavior 
would be to accurately monitor the period history after the 
event. We note that the lack of the recovery of the P value 
to the pre-event one (in contrast with the known behavior of 
glitches observed in radio pulsars) would argue against "mag- 
netar" models, at least in the current formulation, while its de- 
tection might be not conclusive for any model (magnetar and 
accretion). 

1RXS J 170849-4009 10 is also the first AXP for which spec- 
tral changes as a function of pulse phase have been sig- 
nificantly detected. These spectral/timing properties make 
1RXS J 170849-4009 10 an especially interesting AXP to study. 
More sensitive and/or higher throughput observations of 
1RXS J 170849-4009 10 might yield important additional infor- 
mation on the spectral changes causing the pulse shape vari- 
ations and extend the energy range over which the source is 
detected above lOkeV. 
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